Two regimes in the magnetic field response of superconducting MgB 2 
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Using Scanning Tunneling Microscope at low temperature we explore the superconducting phase 
diagram in the 7r-band of the two-band superconductor MgB2. In this band the peculiar shape of 
the local tunneling spectra and their dynamics in the magnetic field reveal the complex character 
of the quasiparticle density of states (DOS). The gap in the DOS is first rapidly filled with states 
in raising the magnetic field up to 0.5 T and then slowly approaches the normal state value : The 
gap is observed up to 2 T. Such a change in the DOS dynamics suggests the existence of two terms 
in the DOS of the 7r-band: a first one, reflecting an intrinsic superconductivity in the band and a 
second one, originating from an inter-band coupling to the er-band. Our findings allow a deeper 
understanding of the unique phase diagram of MgE>2 . 



The discovery in January 2001 of superconductivity at 
39K in the well-known binary compound MgB 2 has 
not only opened a new way for potential applications of 
superconductors but also offered to the scientific com- 
munity an unexpected field of fundamental research. A 
huge effort was very soon rewarded: specific heat 0, 
tunneling 0,0,0 an d point-contact spectroscopy (PCS) 
experiments along with the theoretical predictions 
have shown MgB 2 to be a very unusual supercon- 
ductor in which two electronic bands contribute to the 
superconductivity in a different way. It was established 
that in the superconducting state two gaps open at the 
Fermi level: a leading one, of about —6-7 meV in the 
two-dimensional er-band and a weak one, of —2-3 meV in 
the three-dimensional 7r-band. 

Since then, the effect was extensively studied in a large 
number of experimental and theoretical works (for a re- 
view see 01)- ^ was shown that the superconducting 
gaps are strongly coupled: Both close at the same criti- 
cal temperature of 39K Thus, in MgB 2 superconduc- 
tivity is not described by two independent condensates 
of Cooper pairs but rather by a single condensate of a 
complex nature. 

Various studies revealed a peculiar magnetic field re- 
sponse of MgB 2 . The study of the vortex lattice by Scan- 
ning Tunneling Microscopy/ Spectroscopy (STM/STS) 
evidenced an unusually large size of the 7r-band vor- 
tex cores. This discovery led the authors to suggest that 
gap in the 7r-band is induced by the intrinsically super- 
conducting er band. Heat capacity measurements have 
shown an anisotropy for H|| c and H|| ab [12|. Remark- 
ably, this anisotropy appears only for H> 0.5 Tesla, with 
a large anisotropy developing at higher fields and leading 
finally to H c2 -3 T (H|| c) and -20 T (H|| ab). Neu- 
tron diffraction experiments revealed a rotation of the 
flux lattice by 30° in a field of - 0.6 - 0.8 T, occur- 



ring simultaneously with a sharp drop in the diffraction 
peak intensity 0|. These results were attributed to the 
suppression of the superconductivity in the 7r-band at 
low fields. Experimentally, the magnetic field response 
of the two bands was studied in numerous PCS reports 
[fl R, H3 . Hif . however the 7r-band gap could be evaluated 
with reasonable accuracy only up to 0.5 T. Theoretical 
works regardingthe mixed state in a two-band super- 
conductor 0, 0i E3 have qualitatively explained the 
above findings as a result of the inter-band coupling. Fi- 
nally, though many experiments revealed a modification 
in the superconducting properties at the magnetic field 
of around 0.5-0.8 T, none of the observed changes were 
unambiguously linked to the evolution of the principal 
superconducting parameters such as the superconduct- 
ing gap(s) and/or the inter-band coupling. 

In this Letter we address the question of the magnetic 
field response of the 7r-band by means of tunneling spec- 
troscopy. Using STM/STS we studied c-axis oriented 
surface of MgB 2 single crystals. We found the 7r-band 
DOS to deviate significantly from the BCS shape even 
at T, and to have a non-trivial field dynamics. We 
observed two distinct regimes: A rapid evolution at low 
fields B < 0.4 T and a slow one at high fields B > 0.6 T, 
separated by a crossover region where the DOS remains 
almost unchanged. A clear gap in the 7r-band was still 
observed at 2 T. Such a dynamics suggests the existence 
of two terms in the 7r-band DOS: A first one, reflecting 
the contribution to the superconducting gap from the 
electron-electron interaction via phonons within the ir- 
band, and a second one, originating from an inter-band 
coupling to the er-band. These two terms behave differ- 
ently in the magnetic field, the first term almost vanish- 
ing in the field of 0.6 T. Our experimental findings show 
the need for an additional theoretical effort in the case 
of superconductivity in MgB 2 . 
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Single crystals of Mgl$2 were grow n by a high pres- 
sure method described elsewhere |20j. The experiments 
were carried out on a low temperature Omicron STM 
at Salerno University and on a home-built STM/STS 
in Paris. The tunneling junctions were achieved by ap- 
proaching mechanically cut Pt/Ir tips to the c-axis ori- 
ented surface of the crystals. Due to the peculiar elec- 
tronic structure of the material, c-axis tunnelling allows 
to probe mainly the quasiparticle density of states (DOS) 
in the 7r-band [TjJ, [ljj . In such a geometry the 7r-band 
component of the local DOS was studied in magnetic field 
up to 3T. The superconducting critical temperature was 
determined locally by measuring the evolution of the tun- 
neling conductance spectra as a function of temperature, 
and was found to be 38. 5K. 

The field dependence of the tunneling conductance 
spectra dI/dV(V) in the 7r-band was obtained by fix- 
ing the STM tip in a selected location of the sample and 
by continuously measuring local I — V tunnelling curves 
while sweeping the magnetic field, the field in all mea- 
surements being parallel to the c-axis. The samples were 
cooled in zero magnetic field. At low temperature, the 
field was slowly increased from zero up to the maximum 
value and reduced to zero. Following, the same procedure 
was repeated for the reversed polarity. The field sweep 
rate was ~ 0.01 — 0.05 T/min and the acquisition time for 
a single spectrum was ~10 msec. During such a short ac- 
quisition time the field change of AB ~ 10~ 6 T is, much 
smaller than the field magnitude B itself, AB/B << 1. 
Thus, the field may be considered as constant for every 
single data set. In type II superconductors however, the 
field sweep induces vortex entry into the sample, which 
modifies the spectra. Indeed, while the STM tip remains 
fixed the vortices move and from time to time appear 
under the tip. Statistically some tunneling spectra re- 
flect the DOS in the vortex core, others correspond the 
DOS between the vortices. Averaging over a large (10 2 - 
10 3 ) number of consecutively acquired spectra results in 
a smooth curve representing the quasiparticle DOS spa- 
tially averaged over the vortex unit cell |2l|. In Fig. la we 
show the dynamics of the tunneling conductance spectra 
(averaged over the vortex lattice unit cell) as a function 
of the magnetic field. The spectra are characterized by 
a gap near zero-bias and by only two of the four coher- 
ence peaks at ± 3 meV, consistent with what is expected 
for the tunneling to the 7r-band and in agreement with 
previous experimental data 0, 0, E3> 0] • The main effect 
of the magnetic field on the tunneling DOS in Fig. la is 
the filling with the states inside the gap. The field dy- 
namics of the DOS at the Fermi level is represented by 
the Zero-Bias Conductance (ZBC) in the tunneling spec- 
tra (Fig. lb). The evolution of the ZBC in the magnetic 
field is unusual : At low fields, the ZBC rises rapidly 
and reaches a value of about 80% of the normal state 
ZBC already at 0.4 T. At higher fields however, the dy- 
namics drastically changes: the filling of states becomes 
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FIG. 1: (a) Evolution of the normalized tunneling conduc- 
tance spectra in the magnetic fields up to 2.3 Tesla at T = 6.6 
K. (b) Evolution of the Zero Bias Conductance in increasing 
(triangles) and decreasing (circles) magnetic field. 



much more slower and, even at 2.0T, it is still possible 
to distinguish the quasiparticle peaks and a minimum in 
between in the raw tunneling conductance spectra. At 
intermediate fields (0.4 < B < 0.7 T) there is a crossover 
region in which the ZBC remains roughly constant. Re- 
markably, this range matches exactly the range in which 
a rotation of the vortex lattice was observed [13J . By low- 
ering the field we observe a similar dynamics, though the 
crossover region is slightly shifted to higher fields. Such a 
slightly hysteretic behavior indicates different vortex dy- 
namics in increasing and decreasing fields, which may be 
due to geometrical barriers, vortex pinning, and lattice 
re-arrangements. In order to follow such an unusual field 
dynamics we use two different models describing the DOS 
in the 7r-band. The first one considers the superconduct- 
ing gap to be composed of two terms: a first one origi- 
nating from the electron-electron interaction via phonons 
inside the 7r-band, and a second one, arising from the 
superconducting coupling to the cx-band via phonon ex- 
change. The resulting superconducting gap A x is calcu- 
lated self-consistently ^tJi but finally, the DOS in the 
7r-band remains BCS-like: 
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FIG. 2: (a) Zero field spectrum fitted by BCS model (1): 
T= 7 K, A=1.96 meV and Y Dynea = 0.32 meV. (b) Magnetic 
field dynamics of superconducting gap (black circles) and pair 
breaking parameter V Dynes (open triangles) extracted from 
BCS fit of the tunneling spectra, (c) Zero field spectrum 
fitted using model (2): T= 6.6 K, A° = f.f meV, A° = 6.0 
meV, T n =1.3 meV, r CT = 0.9 meV and Yoynes = 0.3 meV. 
(d) Magnetic field dynamics of superconducting gap (black 
circles) and pair breaking parameter Y Dynes (open triangles) 
extracted from model (2). 



In this Dynes formula |22j the pair-breaking parameter 
r accounts phenomenologically for the effect of the mag- 
netic field on the DOS. 

The second model, in addition, considers the quasipar- 
ticle inter-band scattering. To analyze the shape of the 
DOS in such a case we use, as a first step, two coupled 
equations proposed by McMillan |2^| for the proximity 
effect in real space. This formalism was already success- 
fully applied in some previous works 0> E3- an d it's 
basic idea is well described in |25|. In the framework of 
this approach, the DOS in each band deviates from the 
BCS shape as both A„. and A a become energy dependent 
and the apparent gap in the quasiparticle spectrum of 
the 7r-band is larger than the self-consistently calculated 
superconducting gap, due to the quasiparticle scattering 
from the c-band: 
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The 7r-band DOS is still given by Eq. (1), however is 
replaced by A n (E), whereA° and A° represent the in- 
trinsic pairing potentials. The effect of finite temperature 
is accounted for in both models in a standard way by the 
convolution integral of the DOS (I) with the derivative 
of the Fermi-Dirac function. 

We now focus on the dynamics of the superconduct- 
ing gap in the magnetic field using the above models. 
As we will show in both cases the superconducting gap 
A w exhibits a strong change near 0.6T (for the decreas- 
ing field branch). In Fig. 2a we present an example of a 
BCS fit (1) to the raw tunneling zero field conductance 
spectrum. We note that this fit fails to reproduce simul- 
taneously the value of the ZBC and the amplitude of the 
quasiparticle peaks for any field. However, it does allow 
an estimate for the energy scale of the gap. Fig. 2b shows 
the field dynamics of the main parameters extracted from 
the BCS fits. In particular, one can see that the super- 
conducting gap exhibits a rapid drop at 0.6T while the 
pair-breaking rises continuously in this field region. In 
Fig. 2c we show a typical fit of the same DOS as in Fig. 2a 
using McMillan model (2). Within this model we find 
a much better agreement between the experimental data 
and theoretical curves. Also in this case, the supercon- 
ducting gap extracted from the fits exhibits a remarkably 
strong damping at a field around 0.6T (black circles in 
Fig. 2). In all calculations T w and T a are kept constant 
since physically the inter band scattering rates should 
not change significantly in the magnetic field. 

The use of McMillan model is not arbitrary but mo- 
tivated by the fact that experimentally we observe sys- 
tematic deviations of the tunneling spectra from BCS be- 
haviour. Such discrepancies are better seen in SIS spectra 
as shown in Fig. 3. The spectra were measured at zero 
field using a superconducting MgB 2 tip 4]. In such a 
SIS geometry the spectroscopic features are enhanced by 
the convolution of the DOS of two electrodes. One can 
clearly see a ZB peak characteristic for SIS junctions, 
much stronger quasiparticle peaks appearing at ±6 meV 
and not at ±3 meV as in SIN spectra (Fig.l). Remark- 
ably, some additional features, humps and dips are seen 
on the tails of the quasiparticle peaks. The best fit using 
BCS DOS is plotted as a thin solid line, the thick one 
representing the best fit using McMillan approach (2). It 
is clear that the thick line follows all measured spectro- 
scopic features in finer details. Such an excellent agree- 
ment indicates that one should include into the realistic 
model of the superconductivity in MgB 2 a contribution 
of the quasiparticle inter-band scattering. The values for 
electron-phonon coupling constants should be probably 
reconsidered. 

We see that the use of both models result in a signif- 
icant drop of the superconducting gap energy at 0.6 T. 
This finding is in qualitative agreement with the theoret- 
ical prediction 01 ■ The relative change of the supercon- 
ducting gap is larger within Mac Millan model (2) since 
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crossover region. Our results indicate a rapid suppression 
of the intrinsic term in 7r-band superconductivity for T 
< B < 0.4 T. At high fields (0.7 T < B < 2 T) the super- 
conductivity in the 7r-band survives uniquely due to the 
coupling to the cr-band. The shape of tunneling spectra 
suggests an important role played by the quasiparticle 
inter-band scattering. We think an additional theoreti- 
cal effort is needed to consider properly the inter-band 
coupling in the realistic model of the two-band supercon- 
ductivity in MgB2- 
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FIG. 3: Circles: Zero field S-I-S tunneling spectrum. Thin 
solid line: fit using (1) with A CT = 6.6 meV (relative weight 
2%), A* (98%) = 2.35 meV, V Dynes = 0.1 meV, T = 12 K; 
thick solid line: fit using (2) with A° = 7.4 meV, A° = 1.0 
meV, T CT = 2.3 meV, T n = 2.5 meV, V Dynes = 0.12 meV, T 
= 12 K. Arrows indicate the positions of humps (H) and dips 
(D) perfectly reproduced by the Mac Millan approach. 

there a significant part of the gapped states originates 
from the quasiparticle scattering from the a-band, a pro- 
cess which is field independent. Both models suggest that 
at the field of 0.6T the contribution to the superconduc- 
tivity from the electron-electron interaction via phonons 
in the 7r-band itself is not efficient anymore. At higher 
fields the gap in the tunneling DOS survives due to the 
phonon exchange with cr-band [2(| or due to the both, 
phonon exchange and quasiparticle inter-band scattering 
(Mac Millan model (2) ). 

We note here that while equation (1) is often used to 
describe the quasiparticle spectrum in the magnetic field 
0, its application is not rigorously justified from a theo- 
retical point of view. Indeed, (1) and (2) initially describe 
the quasiparticle spectrum in zero field. Even the intro- 
duction of Tuynes to consider the pair breaking in the 
magnetic field is not really satisfactory as it is correct 
only at B~B c2 . We decided to use such a procedure in 
order to clarify the physical meaning of the field dynam- 
ics we experimentally observed, in the lack of a realistic 
theory describing the DOS in the magnetic field. 

In conclusion, in this paper we studied the supercon- 
ducting 7r-band tunneling DOS in single crystals of MgB 2 
in magnetic field, and we succeeded to observe the su- 
perconducting gap up to 2 T. The evolution of the DOS 
is characterized by two distinct regimes separated by a 
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